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Generalizations of the McMillan map to A/^-body systems 
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Abstract 

The McMillan map is a well-known example of a rational integrable system for one particle in a 
two-dimensional phase space. An elegant recent paper presented a generalization of the McMillan 
map to an A^-body system, for particles moving in d space dimensions. This paper presents some 
alternative generalizations (also completely integrable) of the McMillan map to A^-body systems. 
In all cases, the phase space is foliated by a biquadratic curve in the dynamical variables (and a 
set of suitably chosen angular momentum variables). It is also demonstrated that the constraints 
to generalize the McMillan map to A^-body systems are not trivial. 
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In many cases, the analysis of dynamical systems lends itself naturally to the use of maps 
or so-called Poincare sections. An excellent exposition can be found, e.g., in the text by Dragt 
[1]. (The particle motion around the circumference of a synchrotron provides a natural 
setting for the use of maps.) Of particular interest are maps which describe completely 
integrable systems, both as approximations to real systems and for pedagogical interest in 
dynamical systems theory. The McMillan map [2] is a fully integrable nonlinear system 
in a one-body phase space, and the prototypical example of a so-called rational integrable 
system, i.e. the equations are expressible as rational functions (ratio of two polynomials) of 
the dynamical variables. See [3] and references therein for more details of rational integrable 
systems. An elegant recent paper [4] has presented a generalization of the McMillan map to 
an A^-body system, for particles moving in d space dimensions. 

This paper shows that the generalization in [4] is not unique, and other (also fully inte¬ 
grable) generalizations of the McMillan map to A^-body systems are possible. The present 
paper investigates some alternative models and the conditions for integrability, when gener¬ 
alizing the McMillan map to A^-body systems. The conditions are shown to be not trivial. 
All of the models below are examples of measure-preserving (in fact symplectic) rational 
integrable systems and the phase spaces are foliated by a suitably chosen set of angular 
momenta and biquadratic curves in the dynamical variables, as will be explained below. 

Note that the analysis in [4] was with respect to motion in circular particle accelerators. 
However, the McMillan map is completely general, and the material presented below deals 
exclusively with maps per se, without any specific reference to particle accelerators. 

First some details of the McMillan map [2] are summarized. The notation and terminology 
below mostly follows [3]. For simplicity, it is assumed that all the dynamical variables have 
been scaled to be dimensionless. The “McMillan map” is actually a family of nonlinear 
integrable maps in the (x,p) plane 


f / h't' ' S /-i \ 

X =p, p=-x - . ^ , -. 1 

ap^ + pp + 'f 

Here a, (3, 7 , e and ^ are constants and a prime denotes application of the map. The above 
is an area-preserving rational family of mappings preserving the biquadratic foliation 


ax'^p'^ + f3{x^p -\- xp^) + 7 ( 0 ;^ -I- p^) -\- exp + ^{x -|- p) -f iF = 0 . (2) 


Here iF is a parameter which parametrizes each invariant curve in the plane. Note in passing 
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that the general formalism allows a, /?, etc. above to depend on the invariant K. To avoid 
clumsy notation, this fact will be implicitly assumed below throughout. Following [4], the 
individual particles in an iV-body system will be indexed by a subscript i, hence “Xj” and 
“pj.” An overbar is employed to denote an average, i.e. / = (1/A^) fi quantity 

/. In terms of the above notation, the map in [4] is given by 


Xi= Pi, 


Pi = -Xi - 




ap^ + 7 


( 3 ) 


The authors in [4] have set /3 = ^ = 0 and 7 = 1 . If 7 7 ^ 0, we can scale (a,/ 3 , 7 ,e, 

l,e/ 7 ,^/ 7 ). If 7 = 0, the proof of integrability in [4] still works, but the motion 
may become unstable. It is noted in [4] that the angular momenta Mij = XiPj — xjPi 
(where i ^ j) are all dynamical invariants. A set of A^ — 1 invariants is given as follows, for 
j = 1,..., A^ — 1 (see eq. (7) in [4]) 




( 4 ) 


2 = 1 


A further invariant is given by (see eq. (13) in [4]) 


I = a {xpY + e xp + 7 (x^ + p^). 


( 5 ) 


It was proved in [4] that the set of N invariants in eqs. (4) and (5) are in involution, i.e. their 
mutual Poisson Brackets vanish, hence the system is completely integrable. The model is 
generalized in [4] to treat the motion of N particles in a d-dimensional space. The invariant 
can be expressed in vector notation as I = a (r ■ p)^ -|- e r ■ p -|- 7 (r^ +P^)- I note in passing 
a more succinct derivation of the invariance of the invariant I presented in [4]: 


I' = a (r' • p')^ -|- e r' ■ p' -I- 7 (r'^ -I- p'^) 

-2 — 


= ap- [ r + ^= - 

ap"' -|- 7 


ep 

ep- 1 r + 

ap"' -|- 7 

2 ep^ r ■ p e^p2^ 

= a \ r ■ p~ -\ -=-h — 


ep 


7P^ + 7 r + ^= 


ap^ + 'y 


ap'^ + 'y {ap'^ + '^y 

_ e p^ \ — / 

— e I r ■ p H-=- + 7 p^ -|- 7 -|- 


a p + 7 / 

= J+(2 eFTp+^' 


2 , 2 e r ■ p 


e^ p 2 


ap^ + 7 {ap"^ + 7)" 


ap +7/ Vap^-l-7 


«P^ + 7 _ ^ I 


=0 


( 6 ) 
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Hence I' = I and the invariance is established. 

Consider now, provisionally, a different interaction, also nonlinear, where we square the 

sum of the momenta (as opposed to the sum of squares) 

/ , Ppl + (^Pi + ^ 

Xi=Pi, Pi = -Xi - ^ ^ , — • 7 

a + p p + 7 

This contains too many terms. The angular momenta are not conserved, in general, under 
the application of the above map 

M'j = x[p'j - x'jp'i 


= -pi + 

= Mij + 


Pp] + ePj + ^ 


^ ' ap'^ + [3 p + 'y 
^PiPj - ^ 


+ Pj Xi + 


Ppf + epi + ^ 

ap^ + /dp + 7 


( 8 ) 


{Pi -Pj)- 


ap^ + (3 p + 'y 

The angular momentum is conserved only if /5 = 0 and = 0- Note in passing that a similar 
derivation also explains why the map in [4] also sets /3 = ^ = 0. The McMillan map permits 
nonzero values of jS and ^ because it treats only iV = l, soi = j = l and pi — pj = 0. Hence 
we postulate the following model, where we generalize immediately to a model of N particles 
moving in d space dimensions. With an obvious notation, 

ep 


ri= Pi, 


Pi = -ri - 


+ 7 


(9) 


Here the indices i and j run from 1 to N, and each vector r* and Pi has d components. 
The new definition of the angular momenta is Mijki = {ri)k{pj)i — {rj)i{pi)k. Here k and 
^ run from 1 through d. If i = j then we must have k ytz £ and ii k = £ then we must have 
i j- Note that the above dehnition includes terms such as Xipx 2 — X 2 Pxi, which we do not 
not normally think of as angular momentum. (The same remark applies to the definition 
of the angular momenta in [4].) It is straightforward to derive a set of dN — 1 invariants 
analogous to those in eq. (4), say (where u = 1,..., dN — 1 and we require a sum with 
a messy collection of indices), which will not be written out explicitly here. To obtain the 
hnal invariant, we sum over i to obtain map equations for the centroid {r,p) 

ep 


r =p, 


p = —r — 


+ 7 


( 10 ) 


This is effectively a single particle mapping analogous to eq. (1). It follows from eq. (2) that 
the invariant is 

J = a+ er ■ p + 7 (f^ + p^). (11) 
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Next, the Poisson Bracket of J with any angular momentum vanishes, because J is sym¬ 
metric under an interchange whereas Mijkc is antisymmetric, hence 

{J, Mijki] = 0. This establishes that {J, 1^} = 0 for all z/ = 1,..., dN — 1. Hence the system 
is completely integrable. It is immediate that we can generalize a to a vector: 


ri =Pi 


P'i = 


epi 


( 12 ) 


(X ■ p + J 

The angular momenta Mijki are all conserved, hence also the /j,. The additional invariant is 


j = cxr^ cxp^ + ea.ra.p + ^ {a. r^ + cx - p^). 


(13) 


It is straightforward to derive the proof of invariance. We cannot, however, generalize e to a 
tensor pi —)■ J2m^imPm, because the angular momenta will not be map invariants. For the 
map in [4], we can go one step further and generalize a to a symmetric tensor. With an 
obvious notation, the map equations are 


ri = Pi, 


Pi = -ri - 


^Pi 


(14) 


p ■ cx ■ p + 

The angular momenta Mijki are map invariants, hence also the ly. The additional invariant 
is, written in an explicitly symmetric form 


i = rcxrp(xp + ^e{rcxp + pcxr) + 'y{r(xr+p(xp). 


(16) 


The above are all models of completely integrable systems. Let us consider some other 
related models and investigate the invariants. Consider a map with constant vectors as 
follows: 


ri= Pi, 


Pi 


= -Vi 


epi + ^i 
+ 7 


(16) 


Then 


M-jke = (riY^i - 

{Pj)i + i^j)£ 
+ 7 

{Pi)k{^j)^ {Pj'jii^ijk 


= -(Pi)fc( 

= Mijki + 


+ {Pj)i[ {rYk - 


^ iPi)k T 

ap^ + 7 


(17) 


ap‘^ + 


Hence the angular momenta are not map invariants. Once again we sum over i to obtain 


centroid equations 


r' = p. 


P 


ep + l 

ap^ + 7 


(18) 
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This is again effectively a single particle mapping analogous to eq. (1) and it follows from 
eq. (2) that the invariant is (the invariance is also easy to establish directly) 

J^ = a + er ■ p + 7 (r^ + ■ (h + p). (19) 


Hence for this model is a map invariant but regrettably the angular momenta are not 
conserved. Next, let us introduce a constant vector /3 into the map equations as follows 


r 


/ 

i 


= Pi, 


fBppi + epi 

ap'^ + f3 ■ p + 'y 


( 20 ) 


Then 


M-jke = {riYkiPjYi - {rjYiiPiYk 


= -{Pi)k 



(/3-p + e) (pj)A 

ap"^ + f3 ■ p + 'j J 


+ iPj)i 



(/3 • p + e) {pi)k \ 

Q!p2 + /3-p + 7/ 


^ijki ■ 


( 21 ) 


Hence the angular momenta are map invariants. The most general form of an additional 
invariant, biquadratic in the coordinates and momenta and symmetric under the interchange 
Vi o Pi, is 


JjS = a r^p^ + \{(3 ■ r + (3 ■ p)r ■ p + fi{f3 ■ pr^ + (3 ■ r p^) + e r ■ p + 7 (r^ + p^). (22) 

Here A and /i are constants. Comparing with the McMillan map, we must have A + /i = 1. 
After tedious but straightforward algebra, we obtain 



Jp + 


2(1-A) 


A(/3-p + t) 

a p2 + /3 ■ p + 7 


{(3 ■ pr ■ p — (3 ■ rp^). 


(23) 


There is no value of A which will make the term in the brackets vanish. Next, for iV-body 
motion in one dimension, the vectors reduce to scalars and f3 pf p—f3 f p^ = (3pxp—(3xp'^ = 
0. In 2 or 3 space dimensions, note that/3-pr-p — /3-'rp^ = {px f3) xp)- This vanishes if 
/3|| r xp. Note that r x p is a map invariant because it is a sum of angular momenta (which 
are all individually map invariants) and a sum of invariants is also an invariant. However, 
there is a weak point in that if /3 || r xp, then in the original map equations (3-p = 0, which 
means (3 does not appear in the map equations in the hrst place. Hence an “invariant with 
a term in /3” is vacuous in 2 or 3 space dimensions. In four or more hyperspace dimensions, 
the concept of a vector cross-product is not well-dehned (although the above dot products 
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are well defined). Hence for iV > 1, it is nnclear that an invariant with a term in (3 exists 
for motion in more than one space dimension. As for the model in [4], we can write the map 
equations 


ri = Pi 


Pi = 


f3ppi + epi 

ap^ + (3 ■ p + '^ 


(24) 


The candidate invariant is (again with A + p = 1) 


Ifi = a{r ■ pY + A(/3 ■ r + (3 ■ p) r ■ p + fi{/3 ■ pr'^ + (3 ■ r p^) + e r • p + 7 (r^ + p 2 ). (25) 


Application of the map yields 


= + 


2(1-A)- 


A (/3 ■ p + e) 


{f3 ■ p r ■ p — 13 ■ r p'^ 


(26) 


a p2 + /3 ■ p + 7 

In this case, even for motion in one space dimension, 1^ is not invariant if A^ >1. We must 
have /3 = 0 except for the McMillan map N = d = 1. 

To summarize, the McMillan map [2] is a nontrivial one-dimensional rational integrable 
system. A generalization to A^-body systems (in d space dimensions) was found in [4]. 
This paper presented alternative generalizations of the McMillan map to A^-body systems. 
In all cases, the models are symplectic rational integrable systems and the phase spaces are 
foliated by a suitably chosen set of angular momenta and biquadratic curves in the dynamical 
variables. Additional models with a partial set of dynamical invariants were also displayed, 
i.e. not fully integrable, which demonstrate that there are nontrivial restrictions as to how 
the McMillan map can be extended to A^-body systems in d dimensions. 
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